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The momentum distribution of electrons from semi-leptonic decays of charm and bottom for mid- 
rapidity 1 2/ 1 < 0.35 in p-\-p colhsions at -^s — 200 GeV is measured by the PHENIX experiment at the 
Relativistic Heavy Ion Colhder (RHIC) over the transverse momentum range 2 < px < 7 GeV/c. 
The ratio of the yield of electrons from bottom to that from charm is presented. The ratio is 
determined using partial D/D — > e'^K^X {K unidentified) reconstruction. It is found that the 
yield of electrons from bottom becomes significant above 4 GeV/c in pT- A fixed-order-plus-next- 
to-leading-log (FONLL) perturbative quantum chromodynamics (pQCD) calculation agrees with the 
data within the theoretical and experimental uncertainties. The extracted total bottom production 
cross section at this energy is a^j, — 3.2^ij j (stat)ti'3(sys)/ib. 

PACS numbers: 13.85.Qk, 13.20.Fc, 13.20.He, 25.75.Dw 



Measurements of heavy flavor production (charm and 
bottom) in p-\-p collisions present stringent tests for 
pQCD calculations. For instance, while bottom produc- 
tion at the Tevatron is well described by next-to-leading 
order (NLO) pQCD [H, the cross section for charm pro- 
duction at high pt, though compatible within the the- 
oretical uncertainties, is higher than the preferred the- 
oretical value by ^--^50% 0. Measurement of heavy fla- 
vor in p-\-p collisions also provides an important baseline 
for study of the medium created in relativistic heavy- 
ion collisions. The PHENIX experiment at RHIC has 
measured single electrons from the semi-leptonic decay 
of heavy flavor at mid-rapidity in p-\-p and Au-|-Au col- 
lisions at Y^sjvlv = 200 GeV [3| . Strong suppression 
of the single electron yield at high px, which includes 
contributions from both charm and bottom decays, was 
observed in central Au-|-Au collisions This effect is 
conventionally attributed to energy loss by the parent 
parton in the medium Q ; one also expects the energy loss 
suffered by bottom quarks to be significantly less than 
that suffered by charm quarks due to the difference in 
their masses [1,01. Clearly, for both pQCD comparisons 
and the heavy-ion reference, one wants to disentangle the 
yields of charm and bottom at RHIC energies. 

In this Letter, we present the the yield ratio of single 
electrons from bottom to those from heavy flavor at mid- 
rapidity in p-\-p collisions at ^/s = 200 GeV, using partial 
D/D e^K^X [K unidentified) reconstruction. 

The data were collected with the PHENIX detector [8] 
in the 2005 and 2006 RHIC runs using its two central 
arm spectrometers. Each spectrometer covers \t}\ < 0.35 
in pseudorapidity and Acjj = 7r/2 in azimuth. The arms 



include drift chambers (DC) and pad chambers (PCI, 2, 3) 
for charged particle tracking, a ring imaging Cerenkov de- 
tector (RICH) and an electromagnetic calorimeter (EM- 
Cal) for electron identification and triggering. Beam- 
beam counters (BBCs), covering pseudorapidity 3.1 < 
1 77 1 < 3.9, measure the position of the collision vertex 
along the beam (zvtx) and provide the interaction trig- 
ger. In the 2005 run, helium bags were placed in the 
space between the beam pipe and DC to reduce photon 
conversions. The bags were removed in 2006. 

Two data sets are used for the analysis: (1) a min- 
imum bias (MB) data set recorded with the BBC trig- 
ger, and (2) an electron enriched sample, recorded with a 
level-1 'ERT' trigger requiring a combination of EMCal 
and RICH information in coincidence with the BBC trig- 
ger. The BBC trigger cross section is 23.0 ± 2.2 mb [9|]. 
Since only ^ 53% of inelastic p-\-p collisions satisfy the 
BBC trigger condition, only a fraction of the inclusive 
electron production events are triggered. This fraction, 
which is Pt and process independent, is determined to 
be 0.79 ± 0.02. After selection of good runs and a vertex 
cut of Izvtxl < 25 cm, an integrated luminosity (/ Ldt) in 
the ERT data of 1.77 pb"! in the 2005 run and 4.22 pb"^ 
in the 2006 run are used for this analysis. 

Charged particle tracks are reconstructed using the DC 
and PCI. The momentum resolution is 1% at ~ 
1 GeV/c, and the momentum scale is calibrated within 
1%. Electron identification (elD) is performed using the 
RICH and EMCal. The purity of the electron sample is 
better than 99% for 1 < px < 5 GeV/c 0]. Our previous 
measurement Q determined the spectrum of the single 
electrons from heavy fiavor in the 2005 run. Inclusive 



4 



electron spectra from the 2005 run and the 2006 run are 
consistent within 5% after taking into account a contri- 
bution from the increased photon conversion due to the 
absence of the hehum bags. 

The spectrum of the single electrons from heavy fla- 
vor is determined using the "cocktail method" [1, 0]. 
The electron spectrum from all known sources except 
semi-leptonic decay of heavy flavor is calculated using 
a Monte Carlo simulation and subtracted from the in- 
clusive spectrum in the cocktail method. The dominant 
source of background is the 7r° Dalitz decay. The cock- 
tail also includes contributions from quarkonium {J/tp^ 
T) and the Drell-Yan process, which were neglected 
in our previous measurements Ql- These contribu- 
tions are negligible (smaller than 1% in background) for 
Pt <1 GeV/c but, become significant at high (above 
10% for pT > 2.5 GeV/c) 0. The signal to background 
ratio increases with increasing , approaching unity for 
PT 3 GeV/c 0. 

The systematic uncertainties of the inclusive electron 
spectrum includes the uncertainty in luminosity (9.6%), 
geometrical acceptance (3%), elD efficiency (2%), and 
the ERT trigger efficiency (4% at px > 2 GeV/c). The 
uncertainty in the cocktail method is px dependent (3% 
at px ~ 2 GeV/c, increasing to 9% at 9 GcV/c). 

The ratio of (6 ^ e) to (c ^ e 4- 6 --^ e) is extracted 
from the correlation between the heavy flavor electrons 
and associated hadrons The extraction is based 

on partial reconstruction of the D/D — e^K^X de- 
cay. The invariant mass of unlike charge-sign electron- 
hadron pairs reveals a correlated signal below the D me- 
son mass of ~1.9 GeV/c^, because of the charge cor- 
relation in the D decays. Pairs are formed between a 
trigger electron (2.0<px<7.0 GeV/c) and an oppositely 
charged hadron (0.4<px<5.0 GeV/c). The acceptances 
of positive and negative charged particles are forced to 
be identical by a geometrical acceptance cut. Since the 
momentum range of good charged kaon identification is 
limited, K identification is not performed but the mass of 
all reconstructed hadrons is set to be that of the K. Most 
e^e~ pairs are then removed by an electron veto cut for 
the hadrons. The reconstructed mass of e+e~ pairs has a 
clear peak at low mass. The remaining background e^e~ 
pairs are removed by requiring M^e > 80 MeV/c^, where 
the pair mass is calculated assuming both particles in the 
pair are electrons. 

Depending on the origin of the trigger electrons, the 
inclusive reconstructed electron-hadron pairs are: (1) 
unlike-sign pairs from charm, (2) unlike-sign pairs from 
bottom, (3) combinatorial background where the electron 
is a background electron and (4) background from unlike- 
sign hadron-hadron pairs due to hadron contamination in 
the electrons. The main background source is the combi- 
natorial background (3) and almost all background elec- 
trons are from e~^e~ pair creation. Like-sign electron- 
hadron pairs are used to subtract this background. Since 



electrons from e+e~ pair creation and hadrons do not 
contribute to charge correlated signals, subtraction us- 
ing like-sign pairs cancels out completely the combinato- 
rial background where the trigger electron is from e~^e~ 
pair creation (3). Only the neghgibly small (<1%) con- 
tribution from ifgg decay is not canceled out by the sub- 
traction in the background (3). The contribution from 
hadron contamination (4) is also less than a 1% effect 
due to the excellent electron identification. After the sub- 
traction, the reconstructed pairs include a contribution 
from bottom (2) due to not identifying K. The contri- 
bution from bottom (2) is much smaller than that from 
charm (1) due to the bottom decay modes and kinemat- 
ics. The reconstructed pairs also contain a signal from 
partial reconstruction of heavy fiavor hadrons and a con- 
tribution from a combination of heavy flavor electrons 
and hadrons from jet fragmentation. The ratio of the 
yield of unlike-sign pairs to that of like-sign pairs is about 
1.1 for invariant masses {Msk) below 1.9 GeV/c^. 

The fraction of bottom contribution to the electrons 
from heavy flavor is obtained as follows: 
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where edata is the tagging efficiency in real data and ec{b) 
is the tagging efficiency for charm (bottom) production. 
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tracted unlike-sign electron-hadron pairs for invariant 
mass within 0.4 < Mex < 1-9 GeV/c^, which corre- 
sponds to the mass range of charmed hadrons. Here, 
Nc(b)-ttag is the number of reconstructed signals within 
0.4 < MeK < 1.9 GeV/c^ for charm (bottom) produc- 
tion. 

Figure [T] shows the M^k distribution of the recon- 
structed signals, which is normalized by the yield of 
heavy flavor electrons (A^e(_f/F)) in the range 3 < px < 
4 GeV/c (panel a) and 4<px<5 GeV/c (panel b). The 
tagging efficiency in real data, Cdata, is determined by 
the integration of the M^k distribution in Fig. [1] from 
MeK — 0.4 to 1.9 GeV/c^ as a function of electron px- 

The tagging efficiencies for charm and bottom produc- 
tion, ec and €{,, are calculated with the combination of 
PYTHIA and EVTGEN [H, [l3| . PYTHIA is used to simulate 



charm and bottom production in p+p collisions at ^/s = 
200 GeV and is tuned to reproduce heavy flavor hadron 
ratios: D+/D° = 0.45 ± 0.10, D,/Do = 0.25 ± 0.10, 
Ae/D° = 0.10 ± 0.05, B+/B° = 0.5a_B^% = 0.40 ± 



0.20, and Sbaryon/B° = 0.20 ± 0.15 0, Q, [il, H [Hi 
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FIG. 1: (color online) Comparison of data to a pythia and 
EVTGEN simulation of the invariant mass distributions in 
PHENIX acceptance for the reconstructed signal in the 2006 
run. The electron px range is 3.0 - 4.0 GeV/c (a) and 4.0 - 
5.0 GeV/c (b). The ratios, {b ^ e) / [c ^ e + b ^ e), in solid 
Unes are 0.26 (a) and 0.63 (b). Error bars (boxes) indicate 
statistical (systematic) uncertainties. 



EVTGEN, which is a Montc-Carlo simulation suited for 
decays of D and B hadrons, is used to simulate the semi- 
leptonic decays. The dashed (dotted) lines in Fig.[T]show 
the Mf-K distributions of the reconstructed signal for 
the simulated charm (bottom) production for an electron 
3<pt<4 GeV/c (panel a) and 4<pT<5 GeV/c (panel 
b). Some fluctuations in the simulated curves in Fig. [T] 
come from the limited statistics in the simulation, but 
the statistical uncertainties in the simulation are negligi- 
ble compared to that of the data. €c{b) is determined 
in the same way as Cdata from the M^k distribution 
for charm (bottom) production. Since about 85% of 
the extracted signal comes from partcial reconstruction 
of heavy flavor hadrons, the tagging efficiency is deter- 
mined largely by decay kinematics and ^c(b) can be deter- 
mined with good precision. The dot-dash lines in Fig. [1] 
show the contribution from the combination of an elec- 
tron from charm and hadrons from jet fragmentation for 
charm production. The solid line in Fig. [T] shows the 
sum of the MeK distributions for charm and bottom in 
the simulation with the ratio, (6 — > e)/(c ^ e -I- 6 — > e), 
obtained with Eq. [1] 

Systematic uncertainties are categorized into two parts 
related to (1) Cdata in the real data analysis and (2) Cc and 
tb in the simulation study. The dominant uncertainty in 
Edata is the Uncertainty in the number of heavy flavor elec- 
trons (~10%). Uncertainty in Cdata also includes a back- 
ground subtraction uncertainty (1-10%, pt dependent). 
Category (2) includes the uncertainties in geometrical ac- 
ceptance (3%) and the event generator (~ 8% for charm 
and ~ 9% for bottom). The event generator uncertainty 
is based on uncertainties, which are known in the produc- 
tion ratios of heavy flavor hadrons {D^ ID^ , £)g/£)°,etc.), 
known in the branching ratios [3, [H, 13 , estimated 
in the momentum distribution of heavy flavor hadrons 
and estimated in the pythia parameters. 




FIG. 2: (color online) (6 e)/(c ^ e + 6 ^ e) as a function 
of electron pT compared to a FONLL calculation [l^. The 
points show the experimental result. Vertical arrows are used 
to indicate upper and lower limits. The solid line is a FONLL 
prediction and the dotted lines represent the uncertainty of 
this FONLL prediction. 



Figure [2] shows the resulting bottom fraction, (6 
e) / {c ^ e + b e) as a function of electron com- 
pared to a fixed-order-plus-next-to-leading-log perturba- 
tive QCD calculation (FONLL) 18]. In this figure, the 
points show the measured (6 — > e)/(c — s- e + b — s- e). 
For the bins with electron ranges 2 < < 3 and 
5<pt<7 GeV/c, 90% C.L. and mean values are shown. 
The solid line shows the central value of the FONLL pre- 
diction and the dotted lines show its uncertainty. 

In Fig.[3l the single electron spectra for charm and bot- 
tom are measured from the ratio, (5 e)/(c — > e -I- 6 — > 
e), and the spectrum of the electrons from heavy fiavor 
decays. The top panel shows the resulting single elec- 
tron spectra from charm (triangles) and bottom (squares) 
compared to the FONLL predictions . The measured 
spectrum of single electrons (circles) is also shown for 
reference. The middle (bottom) panel shows the ratio 
of the measured cross sections to the FONLL calculation 
for charm (bottom) production. The shaded area shows 
the uncertainty in the FONLL prediction. The larger 
mass makes this uncertainty smaller in the case of bot- 
tom quarks. These calculations agree with the data for 
bottom production. The same is true for charm within 
the theoretical uncertainty with a ratio of data/FONLL 
of ^^2. A similar tendency was obtained at the Tcva- 
tron [EH. 

The electron spectrum from bottom shown in Fig. [3] 
is integrated from px = 3 to 5 GeV/c and gives 
4.8!ti'6(stat)t} g(sys)nb. This spectrum is then extrap- 
olated to pt = using the shape predicted by pQCD. 
PYTHIA with varying intrinsic /ct (1-5 < /ct < 10 GeV/c) 
and FONLL with varying factorization (fj-p) and renor- 
malization (/x/j) scales (0.5 < IIf,rI -I- p|. < 2) 
are used to evaluate the systematic uncertainty (12%) 
to this extrapolation. The extrapolation results in a 
bottom cross section at mid-rapidity of da^i/dy \y=Q= 
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FIG. 3: (color online) (a): Invariant cross sections of electrons 
from charm and bottom with the FONLL calculation [3. (b) 
and (c): The ratios of data points over the FONLL predic- 
tion as a function of electron pt for charm and bottom. The 
shaded area shows the uncertainty in the FONLL prediction. 



energy loss of bottom quarks. 

In conclusion, the ratio of the yield of electrons from 
bottom to that from charm has been measured in p+p 
collisions at ^/s =200 GeV. The ratio provides the first 
measurement of the spectrum of electrons from bottom 
at RHIC. FONLL calculations [l^l agree with this re- 
sult, which provides an important baseline for the study 
of heavy quark production in the hot and dense matter 
created in Au-|-Au collisions. 
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0.92+°;3^(stat)+°;3g(sys)/zb, using ah e total branch- 
ing ratio of 10 ± 1%, calculated using the heavy flavor 
hadron ratios described above. Using hvqmnr [19| with 
CTEQSm [13] parton distribution functions (PDE's) to 
integrate over rapidity, the total bottom cross section is 
determined to be o-^b = 3.2t};i(stat)ti;3(sys)^b. Var- 
ious PDE's and bottom mass values are used to eval- 
uate the systematic uncertainty (8%) of the rapidity 
extrapolation. This result is consistent with our re- 
sult from the dielectron spectrum, which gave df,^ = 
3.9 ± 2.5(stat)l^(sys)^b FONLL predicts a^i = 

1.87^o g7^b, in agreement with both these experimental 
results. 

The fraction of bottom in heavy flavor electrons is 
found to be larger than 0.33 with 90% confidence level 
at Pt > 5 GeV/c. Furthermore, the assumption of 
no bottom suppression directly leads to a lower limit 
on the nuclear modification factor of single electrons, 
i?AA7 of greater than 0.33 with the same confidence 
level. However, according to our measurements, i?AA is 
~ 0.25 ± 0.05(stat) ± 0.05(sys) at 5 < < 6 GeV/c Q 
in the 0-10% central Au-|-Au collisions. At the same time 
the current level of uncertainty in the measurement pre- 
cludes us from placing significant limits on the possible 
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